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Abstract 

We discuss the role of the geometrical phase in predissociation dynamics of vibrational states 
near a conical intersection of two electronic potential surfaces of a D^h molecule. For quantitative 
description of the predissociation driven by the coupling near a conical intersection, we developed 
a method for calculating lifetimes and positions of vibrational predissociated states (Feshbach 
resonances) for X3 molecule. The method takes into account the two coupled three-body potential 
energy surfaces, which are degenerate at the intersection. As an example, we apply the method to 
obtain lifetimes and positions of resonances of predissociated vibrational levels of the 2^ A' electronic 
state of the H3 molecule. The three-body recombination rate coefficient for the H-l-H+H H2+H 
process is estimated. 
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Non-Born-Oppenheimer interaction near a conical intersection (CI) between two potential 
energy surfaces (PES) of polyatomic molecules plays an important role in dynamics of the 
molecules (see, for example, review 1|]). At the exact point of intersection, the interaction 
becomes infinite and the Born-Oppenheimer adiabatic (BOA) approximation breaks down. 
Therefore, the description of the nuclear motion of the molecule with energies around or 
above the energy of a CI should explicitly take into account the two intersecting PESs 
coupled by the non-adiabatic interaction. 

One of the simplest processes involving the interaction near a CI is the H+H2(f',j') — >■ 
H+H2(f, j) scattering process (and its isotopic variants) that has been extensively studied 
in theory and experiment 2, ^. The related processes are the coUisional dissociation of H2 
and its inverse, the three-body recombination of hydrogen: H+H+H^ H2(f , j)+H, which is 
responsible, for example, for the formation of the first generation of stars [4]. The processes 
are governed by the two lowest molecular potential surfaces I'^A' and 2"^ A' of H3 (see Fig. [1]). 
The lowest I'^A' potential is repulsive and leads to the H2+H dissociation. The dissociation 
limit of the 2"^ A' potential correlates with the H+H+H breakup and has a number of quasi- 
bound vibrational levels, which are predissociated towards to the H2+II dissociation due to 
the coupling near the CI. The predissociation of 2'^ A' vibrational states has been studied by 
Kupperman and collaborators [5| using a time-independent scattering framework, the time- 
delay analysis, and a combination of Jacobi and hyperspherical coordinates. The calculation 
has been done with the two-channel diabatic potential of II3 with non-diagonal diabatic 
couplings obtained directly from first derivatives of ab initio BOA electronic wave functions 
of the two interacting states. In another study by Mahapatra and Koppel 6|], the time- 
dependent approach was employed by propagating a wave packet that is initially placed 
on the 2'^A' PES; the lifetimes are then derived from the autocorrelation function. In Ref. 
0, the autho. *o .ed a d.aba.c ..„,o. of the coupled potent.a,. but the 
diabatization is done differently than in Ref. and does not require explicit calculation 
of first derivatives of the BOA electronic wave functions. The lifetimes obtained in the two 
studies are significantly different. 

In this study, we suggest a general theoretical method to describe the nuclear dynamics 
involving two molecular potentials coupled by a non-adiabatic interaction near the CI. As an 
example, we calculate lifetimes of several 2"^ A' predissociated vibrational levels of H3. The 
method can be used for other small polyatomic molecules where the CI plays an important 
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FIG. 1: (Color online) The two lowest ab-initio potential energy surfaces of H3 shown as functions 
of hyperangles < 9 < tt/4: and < (p < In for a fixed hyper-radius p = 2.5 oq. The projection at 
bottom of the figure corresponds to I'^A' PES. 

role. There are two main ingredients in the proposed method: (1) The diabatization proce- 
dure is made in a way that accounts for the Jahn- Teller coupling between the I'^A' and 2'^ A' 
molecular states and the D^h symmetry of the total vibronic wave function. (2) Nuclear dy- 
namics is described by Smith- Whitten hyperspherical coordinates [7], adiabatic separations 
of hyperangles and hyper-radius along with the slow-variable discretization (SVD) 
and a complex absorbing potential (CAP) to obtain resonance lifetimes. 

Treatment of nuclear dynamics. To treat the vibrational dynamics in hyper- 
spherical coordinates p, 9, 0, we represent the vibrational wave function ip{p, 9, 0) as the 
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expansion ip{p,9,(f)) = X^fc ?/fc(P) ^) 0)cfc in the basis of non-orthogonal basis functions 
Vk = T^j{.P)Va,j{Gi<P)- In this expression, {paj{9,(p) is the hyperspherical adiabatic (HSA) 
state calculated at a fixed hyper-radius pj, 

Hp^^PajiO, 4>) = Ua{pj)^a,j{0, 4>) , (l) 

with the corresponding eigenvalue Ua{pj)- Hp. is the two-channel vibrational Hamiltonian 
with the hyper-radius fixed at pj. If pj changes continuously, energies Ua{pj) and the wave 
functions (Pa,j{9, (p) form the HSA curves Ua{p) and channel functions (Pa,p{0, (p), respectively, 
(see Figs. [2], [3l H]). In the basis i/k, the Schrodinger equation for the total wave function 
iplp, 0, (f)) takes the form of a generalized eigenvalue problem for coefficients = Cja 



Cj'a' 



= (^i I ^i' ) C'ja j'a' Cj^a' ■ (2) 

with overlap matrix elements 

Oja^j'a' = {^a,j{0,(l))\^a'j'{0,(f))), (3) 

that replace the familiar non-adiabatic couplings between the ipa,p{d, 4>) channels. This way 
of representing the non-adiabatic hyperspherical couplings provides an important advantage 
over the familiar method of dealing with the couplings using the first and second derivatives 
of Lpa,p{Oi 0) with respect to hyper-radius. In order to obtain lifetimes of predissociated 
vibrational levels, we place a CAP at a large value of the hyper-radius (for details, see Ref. 
jol). The total three-body rotational angular J momentum is in the present calculation. 

Diabatic basis for the coupled H3 potential. We use the adiabatic ab initio l^A' 
and 2'^ A' PES from Ref. [l^, which will be referred as Vi and V2. The dependence of the 
potentials on the two hyperangles is shown in Fig. [H In our approach, we also use a diabatic 
representation of the interaction potential. However, we derive the non-diagonal diabatic 
coupling elements from the ab initio PESs without calculating them explicitly as derivatives 
of Born-Oppenheimer electronic states. This way of representing the non-adiabatic couplings 
has been very successful in diatomic molecules: If the two ab initio PESs have an isolated 
avoided crossing, the vibrational dynamics is well represented by a 2 x 2 diabatic potential 
with the geometry-independent non-diagonal matrix element equal to the half of the splitting 
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FIG. 2: (Color online) Hyperspherical adiabatic potential curves, obtained from uncoupled l^A' 
(black curves) and 2'^ A' (green curves) PESs of H3. Different dissociation limits for the I'^A' family 
correspond to different v and j. Here, we only show the curves of the Ai irreducible representation. 

between the ab initio PESs at the avoided crossing. The diabatic non-diagonal coupling in 
our model is derived in the following way. 

The two I'^A' and 2"^ A' electronic states become degenerate at the equilateral configuration 
and should be referred to as two components of the E' irreducible representation of the D^h 
symmetry group. It is convenient to use the basis functions |-E'+) and \E^) in the two- 
dimensional E' space [ll|, Although the electronic states I'^A' and 2'^ A' for clumped 
nuclei are classified according to the Cs symmetry group, the vibronic states of H3 should 
be classified according to the D^h group. The most general form of the diabatic potential in 
the basis of \E±) for an arbitrary geometry is 

. I A Ce'f\ 
\Ce-^f A j 

where A, C, and / are real- valued functions of the three hyperspherical coordinates 13|. 
The diagonal elements A are the same because of the degeneracy of E± states. The functions 
A and C transform in the D^h symmetry group according to the Ai representation, and / 
has the following property under the C3 symmetry operator: C3/ = / + 27r/3. The A and 



C functions are uniquely determined from the I'^A' and 2"^ A' PESs: A = {Vi + V2)/2 and 
C = (V2 — Vi)/2. The actual form of function / can be derived near the CI: It is equal to 



the phase of the asymmetric normal mode distortion [12|, ll3|. Although this form of / is 
derived near the CI, we will use it everywhere. This is justified because (1) the transition 
between adiabatic states occur only near the CI, (2) far from the CI the phase factor e^*-^ 
does not play a role as long as the symmetry property mentionned above is satisfied. 

The vibrational wave functions ip obtained from Eq. ([2]) have two components ip^ corre- 
sponding to the two diabatic E'^ basis functions. In the adiabatic basis, corresponding to 
the I'^A' and 2'^ A' electronic states, the two components ipi^2 of the ip function have the form 

^2 = i(^-e*^/' - ^l^+e-'^/^)/V2. (5) 



After applying the C3 symmetry operator three times, the molecule returns back to its 
original position, however the components ^/'i^2 change sign ?/'i^2 ~V'i,2 because / /+27r. 
It is a well-known property of adiabatic states in the presence of CI, which is often referred 
as geometrical or Berry phase effect. Because the adiabatic electronic wave functions |1) 
and 1 2) also change sign, the total vibronic wave function 

^ = + 7/>2|2) = + i^-\E-) (6) 

is unchanged after the identity operator C| is applied. 

Results. After solving the hyperangular part of the three-body Hamiltonian, Eq. ([1]), for 
each BOA PES separately, we obtain two uncoupled 'families' of adiabatic potentials Ua{p), 
which are shown in Fig. [21 The curves belonging to the different families can cross. When 
the coupling is turned on between the two electronic states, the crossings turn into avoided 
crossings, which is demonstrated in Fig. [31 The figure shows the HSA curves obtained by 
solving Eq. ([1]) with the diabatic electronic potential of Eq. (jl]). Far from the (avoided) 
crossings the calculations with coupled and uncoupled electronic states produce almost the 
same HSA states. The coupling changes the HSA states at (avoided) crossings only. It is 
also worth to mention, that if the phase factors e^*'^ in Eq. (jlj) are neglected, the dynamics 
described by such diabatic potential is exactly the same as the dynamics with the uncoupled 
BOA PESs. It is because the transformation diagonalizing the operator V is independent of 
the nuclear coordinates if / = 0. Fig. [H shows the vibrational HSA functions ipa,j for several 
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FIG. 3: Close up look at avoided crossings in HSA curves when the PESs are represented by Eq. 

The figure corresponds to the frame shown in Fig. [2l where HSA curves are calculated from 
the uncoupled ab initio PESs. The horizontal dashed lines show the positions of predissociated 
2"^ A' levels. 

adiabatic states calculated with the coupled potential V of Eq. (jll). 

As Figs. [2] and [3] show, the ground and first excited HSA potentials of the 2"^ A' family 
have minima and can have vibrational levels that are pre-dissociated due to the coupling 
with the I'^A' state. Although each vibrational level has components from all the HSA 
curves shown in Fig. [2] (or Fig. [3]), only one component is dominant. Thus, the vibrational 
levels can be characterized by (1) the dominant component a and by (2) the number vi 
of quanta along the hyper-radius. In addition, each HSA curve in the 2"^ A' family can be 
characterized by the number f 2 of quanta in the hyper-angular space and the number I2 of 
those V2 quanta along the cyclic hyperangular coordinate 0. Therefore, each predissociated 
level can be numbered with the triad {f 1, f2^} similar to the normal mode notations for C^y 
molecules. 

Below the H+H+H dissociation limit there are only two {fiO°} and {fi2°} series of Ai 
predissociated levels. Their positions and lifetimes are given in Table [H Lifetimes for the 
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FIG. 4: (Color online) Wave functions of HSA states ipaj of the Ai irreducible representation as 
functions of the hyperangles 9 and for pj = 2.5 qq and several different a. Each wave function has 
two components, (p-^- and (p- corresponding to the two channels of the potential V. The relationship 
between the hyperangles and the three-body configurations they represent is mapped on Fig. 6 of 
Ref. b 



two series are very different. The reason for the difference is that the avoided crossings 
are significantly wider for the O" curve than for the 2° curve (see Fig. [3]). The table also 



compares the obtained results with two other studies 0, 6| of the 2'^ A' predisso dated levels. 
In Ref. js] the non-diagonal couplings in the diabatic basis are obtained from the accurate 
ab initio non-Born-Oppenheimer couplings between the I'^A' and 2"^ A' states, but not from 
PESs as in this study. The agreement with our study for the lifetime and energy of the 
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Er,T; this work 


Er,T; Ref. \^ 


Er,T; Ref. ^ 


{0,0°} 


-3.85, 13. 


n.a. 


-3.74, ~ 3 


{1,00} 


-3.11, 13. 


n.a. 


-3.01, ~ 3 


{2,00} 


-2.4, 14. 


n.a. 


-2.32, n.a. 


{3,0°} 


-1.8, 14. 


n.a. 


-1.70, n.a. 


{4,0°} 


-1.2, 16. 


-1.19, ~ 15. 


-1.14, n.a. 


{5,0°} 


-0.7, 18. 


-0.42, ~ 17. 


—0.65, n.a. 


{0,2°} 


-0.2, 130. 


n.a. 


-0.22, ~ 4.5 



TABLE I: Positions, Er (in units 10 ^ Eh) and hfetimes, r (in fs) of pre-dissociated 2'^ A' vibrational 
levels. Energies are relative to the H(ls) + H(ls) + H(ls) dissociation. 

{40°} resonance is very good as well as for the lifetime of the {50°} level. The agreement for 
the energy of the {50°} level is not as good probably due to the somewhat special character 
of the {50°} level: its hyper-radial wave function extends to relatively large values ~ Toq 
of hyper-radius (and internuclear distances ~ 5ao). In Ref. [6|], the diabatization procedure 
was based on PESs similarly as it is made in the present study, but with one important 



^ are non-equivalent. It means that the D^h 
6| is broken. The lifetimes obtained in Ref. 



difference: the diabatic electronic states in 
character of the vibronic wave functions in Ref. 

6| are significantly different from the present values and values of Ref. js]. The disagreement 
is attributed to the choice of the diabatization procedure in that does not respect the 
degeneracy of the diabatic electronic wave functions. 

It is worthwhile to note the relatively small number Np = 96 of hyper-radial points pj 
required to obtain converged lifetimes and positions of resonances considering the large num- 
ber of sharp avoided crossings (see Fig. [3]). It is an advantage of using the SVD procedure 
rather than the famihar {ipald"^ /dp^\{pa') and {{pa\d/dp\(pa') non-adiabatic couplings between 
the HSA states: It is not necessary to calculate the derivatives on a fine grid of hyper-radius 
near avoided crossings to describe the couplings locally. Instead, in the SVD approach, 
the overlap matrix elements Oia^i'a' between the adiabatic states in Eq. ([3]) account glob- 
ally for the non-adiabatic hyperspherical couplings. It allows us to reduce significantly the 
computational task. 

The obtained lifetimes and positions of the resonances can be used to estimate the three- 
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body recombination rate coefficient for the H+H+H^H2+H reaction. Using the proba- 
bihty P ~ 1/ {tuj) of the transition from 2"^ A' to I'^A', where u is the frequency of oscillation 
in the 2'^ A' potential, and using the formula for from Ref. [l^, taking sum over total 
angular momentum J up to such Jmax that the lowest 2"^ A' vibrational level moves above the 
H+H+H dissociation, and also taking the sum over all three types of vibrational levels {Ai, 
A2, and E), we obtain value ^3 ~ 2 x lO'^'^cm^/s at 300 K with estimated error of about 
50%. This value is in good agreement with other estimation 4], ^3 ~ 2.2 x 10~^°cm^/s at 
300 K, derived using the detailed balance principle from the rate of the H2+H— s>H+H+H 
processes. Rigorous calculations for k^ combining the developed model for the H3 potential, 
the SVD approach, and the three-body R-matrix approach are under way. 

Concluding, we would like to stress that we (1) suggested a model diabatic two-channel 
molecular potential that represents correctly the interaction between vibronic states near 
a CI. The non-diagonal as well as diagonal matrix elements of the potential are extracted 
from the ab initio adiabatic PESs of the molecule. There is no need to use ab initio non- 
adiabatic couplings between PESs to construct the diabatic potential. (2) We applied the 
model potential for calculation of lifetimes of 2'^ A' pre-dissociated levels in H3. (3) Using 
the obtained resonances, we derived the rate coefficient for the three-body recombination of 
three hydrogen atoms. The developed techniques combining the model diabatic potential 
and the numerical method can also be applied to other systems where a CI is expected to 
play a role. 
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